Evidence for the Coexistence of Anisotropic Superconducting Gap and Nonlocal 
Effects in the Non-magnetic Superconductor LuNi 2 B 2 C 
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A study of the dependence of the heat capacity G p (a) on field angle in LUN12B2C reveals an 
anomalous disorder effect. For pure samples, C p (a) exhibits a fourfold variation as the field H < H C 2 
is rotated in the [001] plane, with minima along < 100 > (a = 0). A slightly disordered sample, 
however, develops anomalous secondary minima along < 110 > for fioH > 1 T, leading to an 8-fold 
pattern at 2 K and 1.5 T. The anomalous pattern is discussed in terms of coexisting superconducting 
gap anisotropy and non-local effects. 



Exotic superconductors, defined as those that follow 
the Uemura relation between the superconducting 
transition temperature T c and the magnetic penetra- 
tion depth A, T c oc A~ 2 , have physical properties that 
differ from conventional (BCS) superconductors. High- 
T c cuprates, bismuthates, Chevrel-phases, organic, and 
heavy-fermion superconductors were suggested to consti- 
tute this class. For most of the members, the super- 
conducting gap function is highly anisotropic or has gap 
zeros on the Fermi surface 2] . The rare-earth nickel boro- 
carbides RNi 2 B 2 C (R=Y, Lu, Tm, Er, Ho, and Dy) re- 
main a challenge because they share many features in 
common with the exotic superconductors but, like el- 
emental BCS superconductors, fall below the Uemura 
trend 0. 

Unlike BCS superconductors, which exhibit exponen- 
tial temperature dependence in density-of-states(DOS)- 
dependent quantities below T c , the borocarbides have 
power-law temperature dependences in their low temper- 
ature specific heat 0, NMR 1/T\ 0, and thermal con- 
ductivity , indicating that electronic excitations persist 
even well below T c . Recently, compelling evidence for the 
presence of nodes along < 100 > directions has been re- 
ported both from field-angle thermal conductivity and 
field-angle heat capacity measurements |8j of YNi 2 B 2 C. 

Another interesting feature of the borocarbides is that 
the flux-line lattice (FLL) undergoes a field-driven tran- 
sition from hexagonal to square with increasing magnetic 
field H || [001] llEl EH. Kogan et al. successfully in- 
corporated nonlocal corrections to the London model and 
Fermi surface anisotropy in these materials 01 to de- 
scribe the transition within a BCS scheme, i.e., isotropic 
s-wave gap with electron-phonon coupling. Nonlocal ef- 
fects were also used to explain, without resorting to any 
exotic order parameter, the modulation of the upper crit- 
ical field H C 2(a) 0] and magnetization M(a) [3| with 
the angle a between the magnetic field and the crystal 



axes. The two seemingly irreconcilable viewpoints have 
only increased the confusion about the nature of the or- 
der parameter of the borocarbides. Recently, Nakai et 
al. considered the coexistence of an anisotropic gap and 
nonlocal effects in the borocarbides ^^| and thereby ex- 
plained the reentrant FLL transformation [l6j | in terms of 
the interplay between the two effects. In this report, we 
present further evidence for the coexistence of the gap 
anisotropy and nonlocal effects. In clean systems, the 
gap anisotropy effects dominate while both effects are 
comparable in less clean samples. 

Single crystals of LuNi 2 B 2 C were grown in a Ni 2 B flux 
as described elsewhere |17| . Samples A and C were not 
annealed while sample N was annealed at T — 1000 C° 
for 100 hours under high vacuum. Typical sample di- 
mensions were 1 x 1 x 0.1 mm 3 . The crystal axes of the 
sample were determined by two independent methods, x- 
ray and upper critical field measurements as a function 
of magnetic field angle |l3|, which were consistent with 
each other. The T c 's, determined by the point where 
the steepest drop occurs in the resistive superconduct- 
ing transition, are 15.5, 15.9, and 16.1 K for sample A, 
C and N respectively (not shown). The resistivity at 
T c is 2.34 and 1.44 /zf2-cm for sample A and sample N, 
corresponding to mean free paths of 144.5 and 234 A re- 
spectively. Assuming that 16.1 K is the transition tem- 
perature for a pure sample, sample A with T c of 15.5 K 
is equivalent to 0.8 % of Co dopin g on the Ni site, i.e., 
Lu(Nii_ :c Co^BzC with x = 0.008 [3. The disorder in 
sample A may be associated with defects which can be 
removed by judicious post growth annealing [lflj . 

Using an ac technique 0, we have measured the low 
temperature specific heat of non-magnetic LuNi 2 B 2 C as 
a function of magnetic field intensity and magnetic field 
angle. Three samples of LuNi 2 B 2 C with different T c 's, la- 
beled A, C, and N, were studied and revealed an anoma- 
lous disorder effect. The heat capacity of the samples 
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FIG. 1: (a) Heat capacity of sample C at 2.5 K for H ||[110]. 
(b)-(d) Heat capacity of sample A at 2.5, 4, and 8 K, respec- 
tively for H || [110]. Dashed lines are H 1/2 fit and arrows 
indicate the deviation point from if 1 / 2 . 



with higher T c 's shows a fourfold pattern as a function 
of magnetic field angle, confirming the result from the 
other non-magnetic borocarbide YM2B2C that the su- 
perconducting gap is highly anisotropic with nodes along 
< 100 > |3| . In contrast, the heat capacity of the dis- 
ordered sample with the lowest T c (sample A) shows a 
dramatic change in the field-angle heat capacity above 
0.8 T at 2 K: the maxima along < 110 > split, giving 
rise to minima separated by 7r/4 (eightfold pattern) and, 
further, C p deviates from a square-root field dependence 
at the same field, labeled H s \, The eightfold pattern 
and the deviation from H 1 / 2 dependence in sample A 
are consistent with the coexistence of an anisotropic gap 
and nonlocal effects. 

Fig. 1(a) shows the magnetic field dependence of the 
heat capacity of sample C at 2.5 K with increasing (cir- 
cles) and decreasing (crosses) magnetic field. The dashed 
line represents the least-square fit of Co + b(H — Hq)^ ', 
where Co is zero-field heat capacity and /j-qHq = 0.1 T 
is a fitting parameter that takes account of the Meissner 
effect. The best fit is obtained when (3 = 0.46, namely 
the Volovik effect for nodal superconductors [2(], , and 
is consistent with previous reports on the non-magnetic 



borocarbides Fig. 1(b), (c), and (d) shows the heat 
capacity of sample A at 2.5 K, 4 K, and 8 K, respectively. 
The dashed line is the square-root field dependence and 
the arrows indicate where the data deviate from the fit. 
The deviation field shows systematic increase with tem- 
perature, i.e., 0.8 T at 2.5 K, 1.8 T at 4 K, and no clear 
deviation at 8 K. Above the deviation field, the data fall 
below the H 1 / 2 line. 

Field-angle heat capacity directly measures the change 
in the DOS with magnetic field direction. The DOS of 
a d-wave superconductor has been shown to exhibit a 
fourfold oscillation with field angle against crystal axes 
H2- At T = K, the DOS has a simple form: 



N(E,H,a)/N ~ £> 4 (1 + T| sin2a|), 



(1) 



where D4 is a Doppler-shift induced coefficient and T 
describes the oscillation amplitude. The field-angle sen- 
sitive Doppler effect, arising from the supercurrent flows 
circulating around the vortices, leads to maxima in the 
DOS when field is along gap maxima and DOS minima 
when field is along nodes. A 3D superconductor has a 
much reduced oscillation amplitude (« 6%) compared to 
that of a 2D system (w 40%) due to contributions from 
the out-of-plane component. We note that a similar effect 
is predicted for an (s+g)-wave superconductor |2^| . 

Fig 2(a) and 2(b) show the low-temperature field- angle 
heat capacity of sample C at 2.5 K and sample A at 2 K, 
respectively. The samples were field-cooled to 2 K (or 
2.5 K) and were rotated within the et6-plane by a com- 
puter controlled stepping motor at increments of 3°. The 
heat capacity with increasing and decreasing field angle 
showed reversible behavior for all measured fields, indi- 
cating that flux pinning effect is negligible in our mea- 
surement. Background contributions (Cbkg) from lat- 
tice vibrations and thermometry were subtracted in the 
usual manner 8] and the remaining field-induced heat 
capacity AC = C'totai — Cbkg was analyzed in terms of 
AC (a) = c(l + r|sin2o;|) for pure samples. At low 
fields, there is a clear fourfold oscillation with minima 
along < 100 > for both samples, indicating that the zeros 
of the gap are located along those directions, consistent 
with those of YN12B2C 0,01 ■ The oscillation amplitude 
r is about 4 % for all three samples (sample N is not 
shown, but is essentially the same as sample C); this is 
consistent with the 3D superconductivity in the borocar- 
bides [H. 

At 1 T, surprisingly, the heat capacity of sample A de- 
velops minima along < 110 >, producing two sets of four- 
fold patterns or 8-fold, an effect not observed in either 
sample C or sample N with higher T c 's. The crossover 
field from the fourfold to the (4 + 4) pattern of sample A 
lies between 0.6 and 1 T, which is also the point where the 
heat capacity of sample A deviates from the square-root 
field dependence (see Fig. 1). With increasing field, the 
splitting in sample A gradually disappears and the field- 
angle heat capacity recovers its fourfold pattern above 
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FIG. 2: (a) Normalized heat capacity of sample C at 2.5 K 
and a 4-fold oscillation (solid line), AC(a)/c — 1 + F\ sin 2a \. 
(b) Field-induced heat capacity of sample A at 2 K, AC = 
Ctotai — Cbkg, and a 4+4 pattern (solid line) - see text. 



4 T. We also measured the field-angle heat capacity of 
sample A at 4 K to check if the anomalous peak splitting 
persists at higher temperatures (not shown). The four- 
fold pattern now persists to 1 T, evolving into two sets of 
fourfold patterns above 2 T. The 4 T data at 4 K has a 
shape similar to the 2 T data at 2 K. The crossover field 
H s \ increases with increasing temperature. 

We focus on the fact that the anomalous 8-fold pat- 
tern occurs only at sample A which has half the elec- 
tronic mean free path of the sample N while the T c is 
slightly decreased. According to the nonlocal theory by 
Kogan et al. |l2j, the hexagonal-to-square FLL transi- 
tion depends on the electronic mean free path / and the 
superconducting coherence length £ of the sample. Gam- 
mel et al. found that a mere 9% of Co doping onto the Ni 
site in Lul221 can make the FLL transition field at least 
20 times higher than that of pure matrix for H || [001] 
[25| . Because the FLL transition field for pure sample 
is expected to be small 12]-possibly below our measure- 
ment range the nonlocal effects would not influence the 
field-angle heat capacity of sample N (or C). In con- 
trast, the disorder in sample A is expected to increase 
the transition to a higher field, i.e. to a field relevant in 
the field- angle heat capacity measurement. 

When the magnetic field is rotated within ab-plane, 
the transition field may differ with different field direc- 
tions because of the different nonlocal range, i.e., £/£. 
The two different transition fields can be characterized 
by H s i and H S 2- As a magnetic field rotates within the 
a&-plane for H s i < H < H S 2, the FLL will experience 
a structural change (or distortion), i.e. hexagonal for 
H || [100] and square for H || [110]. Since the boro- 
carbides have nodes on the Fermi surface, the DOS will 
differ depending on the FLL structure [2lJ ■ The negative 
deviation from the H 1 / 2 above H s \ in the heat capacity 
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FIG. 3: Oscillation amplitude T of sample C (circles) and 
sample N (squares). The dashed line is from the 3D nodal 
quasiparticle theory with vf = 1.5 x 10 7 cm/sec 0- The 
bottom panel contrasts the oscillation amplitude due to the 
nonlocal effect (circles) and gap anisotropy (squares). 



of sample A (see Fig. lb) also indicates that the DOS of 
the hexagonal FLL is larger than that of the square FLL. 
The additional FLL anisotropy in the DOS will modulate 
the gap-anisotropy oscillation and leads to a (4+4)-fold 
pattern. 

We hypothesize that the two effects are independent of 
each other and have a form of two cusped 4-fold oscilla- 
tions in the DOS: 

AC(cv) =pl+j32(l + r|sm2a|)(l + 7|sm2(a-45)|), (2) 

where pi and p2 are fitting parameters. The value T 
represents the oscillation due to gap anisotropy in pure 
samples (see Fig. 3). The nonlocal effects give rise to 
a 45°-shifted 4-fold pattern and are accounted for by 7. 
The solid line in Fig. 2(b) is the least square fit of Eq. (2) 
and represents the data very well. The oscillations due 
to the nonlocal effects (7) and the gap anisotropy (r) 
at 2 K are compared as a function of magnetic field at 
the bottom panel of Fig. 3. The FLL effect 7 increases 
sharply above 0.6 T and decreases gradually to zero at 
4 T, indicating that the low field corresponds to H s i and 
the high field to H s2 . 

Fig. 4 summarizes the H — T phase diagram of the dis- 
ordered sample A. Unlike pure samples, it has additional 
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FIG. 4: H — T phase diagram of sample A with lowest T c . 
Hsi is the FLL transition for H ||[110] and H sX for H || [100]. 
The dotted lines are guide to the eye. Triangles and squares 
were sketched to show the corresponding FLL. 

phase lines in the superconducting state where the field- 
angle heat capacity shows a crossover from fourfold to the 
(4+4)-pattern (H s i) or vice versa (H S 2). The increase in 
the crossover fields with increasing temperature is consis- 
tent with the nonlocal effects jl^ , adding strength to our 
viewpoint that the anomalous 8-fold pattern is due to the 
coexistence of nonlocal effects and gap anisotropy. We 
note, however, that the difference between the two tran- 
sition fields H s i and H S 2 is larger than that expected 
within the nonlocal theory For discussion, we as- 

sume that the anisotropy in the FLL transition field be- 
tween < 001 > and < 100 > is similar to that between 
< 110 > and < 100 > because the H C 2 anisotropy is al- 
most same for the two configurations. The observed ratio 
H S 2/ H s i w 4 is larger than the predicted ratio 2 but 
is smaller than the reported ratio of 10 in YN12B2C [2^. 
The difference between experiments and theory may at- 
test that we need to account for both nonlocality and 
anisotropic gap nature of the borocarbides. 

In summary, we have studied the non-magnetic super- 
conductor LUM2B2C via the dependence of heat capacity 
on magnetic field angle and magnetic field. Unlike pure 
samples, a slightly disordered sample A shows a deviation 
from H 1 / 2 and a (4+4)-fold pattern in C p (a) at above 
1 T. The anomalous properties were explained in terms 
of the coexistence of gap anisotropy and nonlocal effects. 
These experiments resolve the apparently irreconcilably 
different views on the nature of the order parameter in 
the non-magnetic borocarbides. 
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